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ABSTRACT 





This study investigates the effects of different charge density rise times upon the 
magnitude of the Cerenkov radiation produced by a semi-infinite electron beam. The 
magnetic field pulses were generated for different rise times and relationships between 
the maximum magnitudes and their associated rise times were obtained. These were 
compzred with theoretical relationships derived from a power series approximation. The 
generated results were close to those predicted by theory for short rise times at short 


radiai distances. For longer rise times, the departure from theory was caused by the 
magnitude of the decay portion of the magnetic pulse. This effect could be mitigated at 
longer ris2 times by increasing the radial distance from the electron beam. 
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I. INTRODUCTION 


A. BACKGROUND 

Cerenkov radiation is the electromagnetic emission caused by the passage through 
a dielectric medium of a charged particle that is traveling faster than the speed of light 
in that dielectric medium. | 

Imagine an electron moving through a dielectric medium relatively slowly. The 
electric field of the moving electron tends to repulse the electrons of nearby atoms while 
attracting the nuclei, creating distortions in the atoms which results in the creation of 
temporary dipoles. Thus, the medium around the electron becomes polarized. Because 
the polarization field around the electron is completely symmetric, there is no resultant 
field at large distances. [Ref. 1: pp. 3-4] 

But if the electron is moving faster than the speed of light in the dielectric medium, 
the polarization field is no longer completely symmetrical. While the azimuthal plane 
symmetry is preserved, there is a resultant dipole field along the axis of travel which is 
apparent at large distances. This field is created at cach path element along the electron’s 
track. Each field radiates a brief electromagnetic pulse before the atoms realign them- 
selves and return to their normal shape. [Ref. I: p. 4] This behavior can be seen in Fig- 
ure |. 

The time development of Cerenkov radiation was introduced by Buskirk and 
Neighbours for application to electron beams. Single charged particles produce a radi- 
ated power that is proportional to the frequency. Electrons in an accelerator bunch ra- 
diate coherently, an effect that more than offsets a single particle increase in radiated 
power with frequency. [Ref. 2: p. 3750] 

This time structure effect is not observable with present technology in S or L band 
linear accelerators due to their relatively high fundamental frequency. Induction accel- 
erators should produce observable results in air for energies greater than 25 MeV due 
to their longer electron bunch structure. [Ref. 2: p. 3753] 

Lyman conducted a preliminary study investigating the magnetic field radiated from 
a passing charge bunch traveling over a finite path. She showed that the magnetic field 
magnitude versus time plots depended not only on the observer's position but on the 


specific time conditions of each case. [Ref 3: p. 3] 
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Figure 1. A Charged Particle Moving Through a Dielectric Medium: (a) slow 
speed (b) fast speed [Ref. 1: p. 4] 


Price expanded on Lyman’s study by investigating how the magnitude and shape of 
the radiated field depended on the observer location relative to the time boundaries de- 
fined by relationships between the arrival times of different parts of a pulse. [Ref. 4: p. 
3] 


B. OBJECTIVIS 

One of the key assumptions in the Cerenkov radiation time development is that the 
current pulse has a ramp-front, i.e. that the electron pulse density increases from zero 
to a maximum constant value in a linear fashion. The end of the pulse is assumed to be 
a mirror image of the ramp-front but the effects of the pulse end hav2 not been consid- 
ered in these studies. 

The primary objective of this work is to investigate the effects of varying the rise 
time of the electron density, i.e. varying the slope of the ramp-front, upon the magnitude 
of the Cerenkov pulse. 

Secondary considerations are to compate the shape and magnitude of Cerenkov 


pulses along lines parallel to the axis of electron beam travel, evaluate the maximum 





Cerenkov magnitude as a function of radial distance from the axis of electron beam 
travel, and plot the magnetic field as a function of time and radial distance. 

















Il. THEORY 


A. RADIATION FIELDS 

To calculate the Cerenkov magnetic and electric fields, one must first determine the 
potentials from the moving charge distribution. The relationship between the fields and 
potentials are given by Equations 2.1 and 2.2 [Ref. 2: p. 3750], which are 


B=VxA, (2.1) 
=p — GA 
=-Vb-—-*, (2.2) 


The speed of light in free space is c, while the speed of light in the medium is c. The 
electron pulse is assumed to move with velocity v in the positive z direction. The charge 
density p, is assumed to be concentrated along the z axis while the charge is assumed to 
move with no change in shape [Ref. 2: p.3750] such that 


Plt.) = p(z,t)(x)6(), (2.3) 
ple.t) = pole — 01). (2.4) 


The charge density p, has units of charge per volume while p and p, have units of 
charge per lenght. 

Using the assumption of Equation 2.3 and applying it to the usual retarded solutions 
of the wave equations, one obtains the potentials [Ref. 2: p. 3750} 


O(r,!) = 4 | = p(t’ t')dz’, (2.5) 
— y =e o gt ? 
Alt.) =< | a alt’ st')az'. (2.6) 


Note that ¢ is the permittivity of the medium. Also, ’ is the retarded time and 
R=r-r’ such that 














lr—r’| 


verte z : (2.7) 


The assumption of Equation 2.4 can now be inserted into the potential and a new 
variable u(z’) = z' — uf’ can be introduced (Ref. 2: p. 3750] so that 


or) =+ | + poluiae', (2.8) 


A(r,t) = rs | L p,(u)dz’. (2.9) 


The new variable u(z’) can be written out explicitly where x, y, and z represent the 
coordinates of an observer or observation point and z’ is the coordinate of the front edge 
of the electron pulse such that 

1 
ee nas 
v{x + y* + (z—2’)}2 
u(2’) = 2! —vt+ ad i ska é (2.10) 

The velocity v is in the positive z direction and A is proportional to v. Therfore since 
A has a component only in the z direction, B only has components in the x and y di- 
rections and can be calculated from [Equation 2.1 [Ref- 2: p. 3751] such that 





ae fet toll) |e |e . 
nat |e R je+e]4 By LPol tae", (2.11) 
ete | Balt) |e 118 ; 
Byes ral R |e C | R ap Weolu)lae’. (2.12) 


The first integrals in Equations 2.11 and 2.12 have a falloff of R~ at large distances 


for radiation and will be ignored. Implementing the u dependence on x and y (Equation 


2.10 where éu/éy = 3/R and Gu/dx = x/R) into these equations gives Equations 2.13 and 
2.14 [Ref. 2: p. 3751] in the form 














2 
yoo, ' 
Bese |S poi ude’, (2.13) 


2 
By =~ a | pe Pou)" (2.14) 


Note that the new term p,'(u) in Equations 2,13 and 2.14 is the derivative of p,(u) 
with respect to the function u, which was defined earlier. Equations 2.14 and 2.15 can 
be combined into one equation. Also a new function s can be defined as s = (x? + yyt 
[Ref. 2: p. 3751] such that 


2 
B= — | = Po'(u)dz'. (2.15) 


B. TIME DEVELOPMENT 

To evaluate Equation 2.15, one must consider the dependence of u upon z’ as given 
in Equation 2.10. If Equation 2.10 is plotted in the u-z’ plane, the first two terms are a 
straight line with unit slope and an intercept that changes with time. The third term is 
a hyperbola which opens in the + uw direction with asymptotic slopes of +v/c. In the 
Cerenkov case where v >c, the curve is as shown in Figure 2, The curve translates 
downward to smaller values of u as time increases. This is due to the time term in 
Equation 2.10. [Ref. 2: p, 3751] 

Certain assumptions about the charge shape and configuration must now be taken 
into account. The charge is assumed to move with no change in shape. This was reflected 
in Equation 2.4: p(z,t) = p,(z — vt). Since p and p, are also functions of u, p(u) = p,(u). 
The charge densitv is assumed to have a ramp-front profile, i.e. the density increases 
from zero to a maximum constant value in a linear fashion. The time that it takes for 
the charge density to reach this constant plateau is called the rise time {Ref. 4: p. 13]. 

Only changing currents (those with a nonzero p,) will contribute to the magnetic 
fields of Equation 2.15. If the charge density has a linear rise to a constant value, then 


the derivative p,’(u) will be a constant valued square pulse of magnitude p,,’ as shown in 
Figure 3. [Ref. 2: p. 375]] 

















Figure 2. The Function u(z’) ’, 


At large negative times, the u(z’) curve ¢, (Equation 2.10) is completely above the 
nonzero portion of the p,’(u) curve so that the contribution from p,’(u) is zero and thus 
the magnetic field is zero. The u(z') curve moves downward as time increases until the 
B pulse begins when u(z’) is tangent (curve ¢, ) to the upper portion of the of the p,’(u) 
pulse. The magnitude of the magnetic ficld increases as u(z‘) moves downward with in- 
creasing time until the u(z’) curve becomes tangent (curve ,) with the lower edge of the 
p, (tt) pulse. The nonzero part of Equation 2.15 has its largest extent at this time--from 
u(z,) to u(z;). At later times, Equation 2.15 breaks into two regions of the 2’ axis and the 
B pulse decreases with increasing time because the extent of the integral in the two re- 
gions decreases as a result of the upward turn of the u(z’) curve. [Ref. 2: p. 3751] 

If the derivative p,'(u) = p,,’ = constant, then it can be taken outside the integral of 


Equation 2.15. Also the values of n= c,/c and f = v/c, can be inserted so that 


B= 16pm! | et dz’. (2.16) 








Figure 3. The Functions p,(u) and p,’(u) over the Pulse Lenght 


The relationship of the charge density to the function u (Figure 4) must be deter- 
mined in order to evaluate p,’. This relationship is 


p 
pa- out py (2.17) 


Note that Au = u, — u, (see Figure 1). 


Taking the derivative of Equation 2.17 with respect to u, one obtains p,,’ such that 


' ’ p 
Po (Ut) = Pm =~ Aa ' (2.18) 





Equation 2.18 is now substituted into Equation 2.16 to give 


p : 
Baap fe | ae, (2.19) 








Figure 4. p as a Function of uw at One End of the Pulse 


Since R? = s? + (z — z')? and Au=va where a is the rise time, Equation 2.19 can be 
integrated over the limits of z, to z» Introducing the function w =z — 2’ and integrating 


gives the result \ 


p _) Ww _), W 
B = —nf? 52 4 tanh) tan}. (2.20) 


Note that w, =z —2/ and w,=z-—z/. Thus, the problem is to find the values of w, and 
wy tRef. 3: p. 14] 


C. LIMITS OF INTEGRATION 

To find the limits of integration z/ and z/ , one must consider two cases. The first 
case is where the u(z’) curve minimum is in the non-zero portion of the p,’(u) function, 
ie. between curves 4, and 4, of Figure 2. Ihe second case is where the u(2’) curve mini- 
mum has passed through the non-zero portion of the g,'(u) function, i.e. lower than the 


curve ¢, of Figure 2. 





The first case produces a single integral where the limits of integration are found by 
solving Equation 2.10 for z’. The function u(z’) = u, where u, is the value of the u curve 
minimum at the upper edge of the p,’(u) function. The result is 

{Ay — B72} + / Be - Ay + °(1 — 8) 
ae Sa (2.21) 
(=f) 
Note that A, = u, + vt. 

The second case produces two integrals where the limits of integration are also 
found by solving Equation 2.10 for z’. Two of the limits come from Equation 2.21 while 
the other two come from where the function u(z’) =u, U, is the value of the u curve 
minimum at the lower edge of the p,’(u) function. The result is Equation 2.21 and 


{Ay ~ B°2} + / B{(2 — A,) + 9°11 — 8} 


zy" = = 5 
(1 ~ p") 


(2.22) 
Note that A, = u, + vt. 

Another way to describe the limits is that they are the z’ values of the intercepts of 
the u(z’) curve with the upper and lower edges of the non-zero part of the ,’'(u) function. 
This is shown in Figure §. 

For the two integrals of the second case, the ordering of the limits is very important. 
The lower and upper limits of the first integral are z,/( min) and z,’( min) respectively 
while the lower and upper limits of the second integral are z,’( max) and z,'( max) re- 
spectively. 


D. RISE TIME DEPENDENCE 

To derive an approximation for the maximum magnetic field dependence upon rise 
time, Equation 2.10 must be developed in a power series. Denoting z,’ as the value of 
z' at which Equation 2.10 has zero slope, the values of z,’, u(z,') , and the second deriv- 
ative of u with respect to z’ [Ref. 2: p. 3751] are 


1 
2 -——— 
sna a1] if (2.23) 


2 _— 
u(z,)=Z—vt+ | => = i - (2.24) 
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Figure 5. The z’ Limits of Integration 
au 1 ¢? } y’ ? 
2H) D1)” = 2D. (2.25) 


D is a constant defined for conveinence of use Thus u can now be expressed as a Taylor 
series about the minimum such that 


u = u(2,’) + D(z’ —2,'). (2.26) 


The limits z,’ and z,’ can be written in terms of the minimum value as z,’ =z,’ + Az’ 
and z,' =z,’ — Az’ where Az’ is the value of z’ ~ z,’. The difference u(z’) — u(z,') = a where 
a is the rise time (or width) of the current derivative pulse p,’ (a can be thought of as a 
distance by multiplying it by v). [Ref. 2: p. 3751] Therefore 


1 
Az! = {-£ ha, (2.27) 


Using this value, the fact that Au = va, and the assumption that s and R are slowly 


varying, Equation 2.19 can be evaluated to give 








p + 
Bnax = 08 iar}? (2.28) 


Thus the maximum magnetic field has a dependence upon rise time to the negative 
one-half power. Therefore, as the rise time increases, the maximum magnetic field de- 
creases, [Ref. 2: p. 3752] 

Two points should be emphasized. Equation 2.20 gives the exact maximum magnetic 
field amplitude for « evaluated at time 4, (see Figure 1), Equation 2.28 gives an approx- 
imation for the maximum magnetic field amplitude as a function of the rise time, giving 
an Opportunity to compare and contrast actual behavior. 


12 














II, CALCULATIONS AND ANALYSIS 


The previous work had identified the exact solution for the maximum magnetic field 
amplitude and an approximation for the maximum as a function of the rise time. Until 
this time, no comparison had been conducted between the two to validate the approxi- 
mation. 

Plots of the magnetic field magnitude versus time were generated from evaluating 
Equation 2.20, using the programs listed in Appendix A. These plots have 200 points 
each. The time increment between points was adjusted so that the maximum magnetic 
field magnitude occurred at the fortieth point of the plot. The increment could be so 
adjusted by knowing that the maximum magnitude would occur at the end of the rise 
time a. The variables of interest in the programs were the rise time a, the radial distance 
s, and the axial distance z, Program runs were made for s = 1.0m, 5.0m, and 10.0m; z 
= 0.0m, 1.0m, and 5.0m; and a = 10.0ps, 20.0ps, and 30.0ps. 

Three general trends were observed from these runs. As a izicreases with s and z 
fixed, the overall magnetic field magnitude decreases. This ig due to the a relation in 
Equation 2.20. As s increases with a and z fixed, the overall magnetic field magnitude 
decreases. This is due to the s~' relation in the arctangent function of Equation 2.20, And 
finally, as z increases with a and s fixed, the overall magnetic field magnitude and shape 
were unchanged except that the pulse was delayed in time more and more as z increased. 
This is because z is merely a function of time, as seen in Equation 2.10. This time delay 
behavior can be seen on the horizontal axes of Figure 6, Figure 7, and Figure 8. 

To further investigate B,,,, plots of the maximum magnetic field magnitude versus 
the rise times associated with each magnitude were generated. This was accomplished 
by the programs listed in Appendix B. The magnetic field as a function cf time was cal- 
culated for 200 points by evaluating Equation 2.20, just as did the programs of Appendix 
A. This was done for a fixed s and a fixed z (z was chosen to be zero for all following 
calculations since it has no effect on the magnitude and simplifies the calculation proc- 
ess). The initial rise time was selected to be 10.0ps. The magnetic field magnitudes were 
run through a DO loop to identify the maximum. The maximum, along with the corre- 
sponding rise time, was stored in an output data file. The rise time was then incremented 


by 5.0ps and the entire procedure is repeated three hundred times. The output data file 
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CERENKOV PULSE (GAUSS) 
0.02 0.04 0.06 0.08 0.10 0.12 0.14 


0.30 


Figure 6. 


B-FIELD VS TIME 


5.0 7.0 9.0 
TIME (SECONDS) 
B yersus t at s = 1.0m, z = 0.0m, and a = 10.0ps 











B=F ELD VS: DIME 
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CERENKOV PULSE (GAUSS) 


0.02 


0.00 


339.00 341.00 343.00 345.00 
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Figure 7. B versus ¢ at s = 1.0m, z = 1.0m, and a = 10.0ps 








B=F tELU VS TIME 


CERENKOV PULSE (GAUSS) 
0.04 0.06 0.08 0.10 0.12 0.14 


0.02 


0.90 


1673. 000 1675.000 1677 .000 1679.00 
TIME (SECONDS) *10 


Figure 8. B versus t at s = 1.0m, z = 5.0m, and a = 10,0ps 











of 300 points is potted by a DISSPLA program incorporated into Appendix B. An ex- 
ample plet can be seen in Figure 9. 

he theoretical approximation of Equation 2.28 predicts that the relationship be- 
twcen the maximum magnetic field magnitude B,,, and the rise time a should be such 
that Bu, oc a5, The best way to examine this beliavior is to apply a log-log scale to the 
output data file used to create Figure 9. If the log-log graph produces a linear relation- 
ship, the slope of the graph will correspond to the exponent of a, The log-log plot 
produces a liaear relationship close to that expected of Equation 2.28 for low and me- 
dium rise times, but the slope begins to decrease as the rise time reaches larger values. 
This behavior can be seen in Figure 10 and an explanation follows. 

The exact solution for the maximum magnetic field magnitude can be calculated by 
evaluating Equation 2.20 where the limits z' (Equation 2.21) are for the time 4, (see 
Figure 1), Since 4.=2,+ a, 2’ becomes a function of a. With z arbitrarily chosen to be 
zero, w/s becomes Ka/s where K is some constant. 

The next step is to examine the arctangent function. For large values of Ka/s, the 
arctangent in Equation 2.20 is approximately constant and large changes in a have little 
effect on the arctangent value. Therefore the arctangent behaves as if there is no a de- 
pendence and B,,,, should be proportional to a7! . But the slope of Figure 10 is tending 
towards zero for high @ values instead of towards -1, the exact opposite of what one 
might think would happen. 

There is one main reason why this behavior occurs. It goes back to the two cases 
for the limits of integration of Equation 2.20. The maximum magnetic field amplitude 
was calculated for the single integral case where the limits of integration were evaluated 
at time #, (see Figure 2). The magnetic field amplitude of the “tail” of the pulse was cal- 
culated for the double integral case where the limits of integration were evaluated for 
times ¢ > 4, (see Figure 5). 

Since these calculations are for the radiative ficlds, the magnitude of the magnetic 
field does not decay to zero but instead eventually decays to some constant value. As the 
rise time increases for a fixed observation point, the magnitude of the “tail” increases 
relative to the maximum magnitude of the pulse. Another way to look at the process is 
that large rise time pulses decay slowly compared to short rise time pulses. For extremely 
large rise times, the maximum amplitude is virtually indistinguishable from the “tail”. 
Note that the maximum amplitude for a short rise time is larger than the maximum rise 
time for a long rise time (irom theory, B,,,, oc a-°*). Therefore as the rise time increases, 
the maximum magnetic ficld amplitude cannot continue to decrease linearly because the 
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Figure 9. B,,,, versus a at s = 1.0m 
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Figure 10. Log B,,, versus Log a at s = 1.0m 











maximum amplitude is approaching a constant value as the rise time increases. This 
behavior is shown in Figure 11 and Figure 12, 

But as s increases, it decreases the value of Kal/s and the arctangent leaves the almost 
constant value range and tends towards the linear response value region. As a result, the 
relationship between B,,,, and a behaves linearly over a larger range of a as s increases. 
This behavior can be seen in Figure 13. 

The data file used in the plots of Figure 9 and Figure 10 was generated at a fixed s 
value of 1.0m. To further investigate this behavior, additional data files needed to be 
calculated and plotted for different values of s to see how much of an effect the radial 
distance might have on the relationship between B,,,, and a. The results are shown in 
Table 1. 


Table 1. B,.. VERSUS a SLOPE VALUES 
slope = -0.341 

slope = -0,452 
slope = -0.467 
slope = -0.475 
























The cancelling of the B,,, oc a7! effect at high a values by the magnitude of the 
m.gnetic pulse decay tail is shown in Figure 14. The effect of reducing the argument of 
the arctangent function by increasing the radial distance s and thus returning the 
arctangent function more to its linear response range is shown in Figure 15. 

The relationship between the maximum magnetic fie'” magnitude B,,,, and the radial 
distance s was graphed. The magnetic field was evaluated using Equation 2.20 and was 
plotted as a function of time for 200 points at a constant s and constant a. The maxi- 
mum magnetic pulse magnitude is determined and stored with its corresponding s value. 
The value of s is incremented and the entire procedure repeated until a 300 point data 
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Figure 11. Slow Decay for a= 1000.0ps at s = 10.0m 
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Figure 12. Fast Decay for a= 100.0ps at s = 10.0m 
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Figure 13. Log B,,,, versus Log a at s = 30.0m 
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Figure 14. B versus ¢ at s = 1.0m and a = 1000.0ps 
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Figure 15, B versus t at s = 30.0m and a = 1000.0ps 
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file has been constructed. This data file is then plotted using DISSPLA. The programs 
for these actions are located in Appendix C. An example plot can be seen in Figure 16. 

The last step is to calculate and plot the three-dimensional relation between time, 
the magnetic field, and radial distance. The magnetic field is calculated versus time for 
50 points using Equation 2.20. The magnitude and corresponding time and radial posi- 
tion are stored in output files. The radial distance is incremented and the entire proce- 
dure is repeated 50 times. The result is three output files containing the coordinates for 
a 50x50 array. A 3-D graphing routine reads these array points, converts them to a 
surface, and plots the surface. The programs for this procedure are listed in Appendix 
D. 

Due to the fact that DISSPLA will only handle square arrays for surface plotting, 
the number of points along the B—1r axis had to be limited to 50 to hold computer 
memory requirements within acceptable limits. This somewhat obscures the detail of the 
3-D plot. Therefore a limited 2-D projection of B~ 1s curves for a few s values was cre- 
ated in Figure 17 to help understand the workings of Figure 18. 

An example plot is displayed in Figure 18. If the graph projection of the B— s plane 
is rotated about the ¢ axis, a representation of the Cerenkov radiation cone in space is 
produced, 

Several values were used in the programs listed in the appendices. Electron energy 
was set at 30 MeV. It must be approximately 25 MeV or higher in order for Cerenkov 
radiation to occur. Current was assumed to be 510 amperes. Any other desired values 
can be used for these constants. Except for the purposes of Figure 7 and Figure 8, the 
observer axial distance z was sct to zero for simplification of calculations. The observer 
radial distance s can never be set to zero except for initialization of incrementation use 
in a DO loop. The function u, was set to zero also for simplification. It has the result 
of causing the function u, to simplify to the expression u, = — va. where a is the rise time. 
Finally, the magnetic field calculations were conducted in MKS units and converted to 
CGS units from the relation: I Tesla ~ 10,000 Gauss. 
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Figure 16.  B,,,, versus s at a= 10.0ps 3 
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Figure 17. B versus ¢ for s = 1.0m, 2.0m, and 3.0m 
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B-FIELD VS TIME VS RADIAL DISTANCE 
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Figure 18. B versus ¢ versus s 














IV. CONCLUSIONS AND RECOMMENDATIONS 


The primary objective was to determine the effect of varying the rise time of the 
electron beam charge density upon the Cerenkov pulse magnitude. The theoretical ap- 
proximation of Equation 2.28 predicts that this behavior should be B,,,, oc a-®5. Actual 
results show that this behavior approaches the theoretical value at large radial distances 
due to the decrease of the argument of the arctangent function. 

This study was conducted for the case of a semi-infinite path lenght. Evaluating this 
behavior for the finite case will have some effect on the expected outcome. 

The magnetic field magnitude falls off as the radial distance increases. This was ex- 
pected from the relationship in Equation 2.20. Since the magnetic field must also travel 
farther as s increases, the pulse is time-delayed relative to a position closer to the 
electron beam. By rotating the 3-D graph of Figure 10 about the time axis, an accurate 
spatial description of the magnetic field pulse is obtained. This figure is the expected 
Cerenkov radiation cone, somewhat analogous to the wake created bv a boat moving 
faster than the current speed of the water. 

The shape and magnitude of the Cerenkov pulse is not affected in any manner by 
changes in the observer's axial position. It is only affected by changes in the observer's 
radial position. This assumes that all other possible variables (energy, rise time, current, 
etc.) are held constant. The tail of such a pulse decays to some constant value which is 
not zero. This effect is because the calculations are for the radiative fields. Calculations 
should be conducted with the non-radiative fields as well. 
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APPENDIX A. CERENKOV PULSE PROGRAM 
PROGRAM PULSEI 


Yaininieteinicinist ininintaicininiaininineiaicinicideisteteiriesiel tel ect iie ieee die ioe iota ieeien: 
* THIS PROGRAM IS LOCATED ON THE NPS MAIN FRAME COMPUTER AND IS * 
* WRITTEN USING WATFOR 77. AFTER ENTERING CERTAIN PARAMETERS, THE * 
* SHAPE OF THE CERENKOV PULSE IS GENERATED. THE VALUES OF THE * 
* PULSE MAGNITUDE POINTS FOR THE GIVEN PARAMETERS ARE STORED IN . 


%* AN OUTPUT FILE FOR LATER USE IN A PLOTTING ROUTINE. 
HINARI TRAIN IAT IAT ITIL ITA TRIN THT IRIAN IIASA IR AAI IIR 


REAL N, U1, U2, BETA, CO, ROE, Al, A2, BPRME 
REAL RISETM, D, DD, El, E2, 211, 21F, 221, Z2F 
REAL Bl, B2, W1, W2, YY, XX, Tl, T2, WPI, WPF 
REAL ZPMIN, SSQ, BPSQ, TINC, IAMP, ENER, ENERJ 
REAL S, Z, C, MO, V, NEG, TPRME, B 

INTEGER I 

DIMENSION TPRME(200), B(200) 

DATA TPRME/200*0.0/, B/200%0. 0/ 


OPEN (UNIT=4 ,FILE='TERMINAL' ) 
OPEN (UNIT=7,FILE='PULPLT DATA A') 


FetedederertetededededeseteseveleleteleLeseledreLeleDeR VTL TANTS AAT RATT ISDN RIND ICANN IIA INDIE 
* N IS THE INDEX OF REFRACTION. IAMP IS THE ACCELERATOR CURRENT * 
% IN AMPS. ENER IS THE ELECTRON ENERGY IN MEV. RISETM IS THE ¥ 
PULSE RISETINE IN SECONDS. S IS THE OBSERVER RADIAL DISTANCE IN * 
METERS. Z IS THE OBSERVER AXIAL DISTANCE IN METERS. CO IS THE * 
VACUUN SPEED OF LIGHT IN METERS PER SECONDS. MO IS THE ELECTRON * 
MASS IN KILOGRAMS. ROE IS THE LINEAR CHARGE DENSITY IN COULOMBS * 
* PER METER. we 
HUTT TTT ITNT LTT RIT TLL LLL DIL LINSEED ELE SDL DDE AD ADA ADAM 


44% 4 


N = 1.000293 

ITAMP = 510.0 
ENER = 30.0 
RISETM = 10. 0OE-12 
§ = 1.0 

Z= 5.0 


CO = 2.99792458E08 
MO = 9.109534E-31 

C = CO/N 

ENERJ = ENER*1, 6021892E-13 

V = CO*SQRTC1. O-( CMO**2. 0) *( COS, 0) /CENERJ**2. 0) )) 
IF (V.LE.C) GO TO 60 

BETA = V/CO 

ROE = IAMP/V 

BPRME = N*BETA 


Seteketeveretededede dete tevere sete retested WNC Re Te Te ICT TTI TCT IIT TT THN TEN I TENT IOE 
% SETTING Ul ARBITRARILY TO ZERO FIXES U2. THIS IS TO SIMPLIFY ve 
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* CALCULATIONS. * 
Heinkivinksciekset tert ekea RIT ITT IA RTT RET IT ITT ITT TR TR IT RTE HITE 


Ul 
U2 


0.0 
-V*RISETM 
FIAT IITIRTIT ATI TAIT ITI IT HITT TIMI IISA TIT IIIA IARI 


* WITH 200 POINTS, TINC WAS SELECTED SO AS FOR THE PULSE MAXIMUM * 
* TO OCCUR AT APPROXIMATELY 20% OF THE GRAPH LENGHT. * 
FRIAR AIR IRIRRTRI AT TRINIHITNITITTRTTITRIERTRTN RTI 


TINC = RISETM/40. 0 


SSQ = S**2.0 

BPSQ = BPRME**2. 0 
D = Z*BPSQ 

DD = 1 - BPSQ 


NEG = SSQ*DD 


Hedereredereredevere dese rere tedede Fee vera kewewe re vewewevewedvewewzewe KVRAF ER ME TE TEN TENE TET 


* ZPMIN IS THE Z VALUE OF THE MINIMUM OF THE U CURVE. * 
deaeaaeinvedesieseden sears ae ar a IT HIT RTE TEN TET INET TEIN TER TEN TE TTT TEM TET TET NIE TE 


ZPMIN = Z = (S/(SQRT(BPSQ - 1.0))) 


Sedvterereverevededededeseye dese de seven reve revere Nese de te Fekete Te Neve Ve HRT NE TE HET TE TENE TE ITE TE NE TET TE TE TE TEE NIE TE TE TET 


%* T1 AND T2 ARE THE TIMES THAT THE U CURVE MINIMUM IS TANGENT TO * 
%* THE Ul AND U2 LEVELS RESPECTIVELY. we 
Severededevededededededetedede sede tevese dese dete TeV NN PLLC TTA IANA ILI IAIN ANSI INNING 


Tl 
T2 


(ZPMIN - U1)/V + (SQRT(SSQ + (Z - ZPMIN)**2.0))/C 
T1 + RISETM 


DO 40 I = 1,200 
IF (1.EQ.1) THEN 


TPRME(I) = Ti 

B(I) = 0.0 
ELSE 

TPRME(I) = Tl + (REAL(I)*TINC) 
END IF 


Seiededededededesedesdevededeacsedesevedetevevededese sede de seve dete se Fee Keke He TTR ITN TIENT IR TLE TTR TT TONITE 


% THIS IF STATEMENT HANDLES THE CHANGEOVER IN LIMITS OF INTEGRAT- * 
ve 


%* JON FROM ONE INTERVAL TO TWO INTERVALS. 
Setetevetedededetedededetesete de Tete LeVTeLelh TL RARER IARI ISELIN SANNA LISSA INS IIA ALN 


IF (TPRME(1).GT.T2) GO TO 10 
Al = Ul + V*TPRME(I) 
= (Z = A1)**2.0 


Pededeteledevete te Ke Pelee Veda re let dK AAT LRT BIRR RESELLER LDA NSAI III IIR. 


* THIS IF STATEMENT HANDLES THE CASE OF HAVING A NEGATIVE SQUARE * 


% ROOT IN REAL SPACE WHEN DETERMINING THE LIMITS OF INTEGRATION. * 
dacdevovedededevererededeve dete sete Tee reve se ge He Fe He I TEN TCT NET TENN TENE TTT TE TE TEETER TEE TE TENE TENET TENE TEINS 


IF (E1. LE. ABS(NEG)) THEN 
B(I) = 0.0 
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GO TO 
are 


0 


3 
((A1-D} + SQRT(BPSQ*(E1+NEG) ))/DD 
((A1-D) - SQRT(BPSQ*(E1+NEG)))/DD 


ENDIF. 


WPI = Z - Zil 

WPF = Z - Z1F 

W1 = WPI/S 

W2 = WPF/S 

YY = ATAN(W1) 

XX = ATAN(W2) 

BCI) = -ROE*N*(BETA**2, 0)*(YY-XX)*1. 0E04/U2 
GO TO 20 


CONTTNUE 

Al ='U1 + V*TPRME(T) 

A2 = U2 + V*TPRME(I) 

El = (Z - Al)**2.0 

E2 = (Z - A2)**2,0 

IF (E1. LE. ABS(NEG)) THEN 
B(I) = 


ul 
NaN 


(A1-D) + SQRT(BPSQ*(E1+NEG)))/DD 
(Al-D) + SQRT(BPSQ*(E1+NEG) ))/DD 


BCI) = 0.0 


((A2-D) + SQRT(BPSQ*(E2+NEG) ))/DD 
((A2-D) - SQRT(BPSQ*(E2+NEG) ))/DD 


WPF/S 

ATAN(W1) 

ATAN(W2) 

~ROE*N*( BETA**2. O)*CYY + XX)%1. OE04/U2 


~ 
wm 
wuktunu 


WPI = Z - Z2F 
=2Z + ZIF 
WPI/S 
WPF/S 
ATAN(W1) 
ATAN(W2) 
B2 = -ROE*N*(BETA**2.0)*(YY - XX)*1. 0E04/U2 
B(I) = Bl + B2 


mx RS 
ws 


oui nn 


CONTINUE 








30 CONTINUE 
40 CONTINUE 


DO 50 I = 1,200 
WRITE(7,1000) TPRME(T) ,B(T) 
1000 FORMAT(G15. 7,615. 7) 
50 CONTINUE 


GO TO 70 

60 WRITE(4, 1010) 

1010 FORMAT(/,'V IS LESS THAN C, THEREFORE THERE IS NO CERENKOV RADIATI 

SON FOR THIS ELECTRON' ,4X,' ENERGY') 

70 CONTINUE 

Cc 

STOP 
END 


PROGRAM PULPLT 


Qa aang 


REAL X(300) ,¥(300) ,XMIN, XMAX, YMIN , YMAX 
INTEGER N 
CHARACTER*20 TITLES 


TITLES = 'B-FIELD VS TIME$' 
N = 200 


READ(7,8000) ((X(I),Y(1)),1=1,N) 
8000 FORMAT(2G15. 7) 


CALL TEK618 

CALL SHERPA('P6 B43) 

CALL PAGE(6,7.5) 

CALL AREA2D(4.5,5.5) 

CALL XNAME('TIME (SECONDS)$' , 100) 

CALL YNAME('CERENKOV PULSE (GAUSS)$' ,100) 
CALL HEADIN(TITLES,100,1.5,1) 

CALL CROSS 

CALL RANGE(X, XMIN, XMAX,N) 

CALL RANGE1(Y,YMIN, YMAX,N) 

CALL GRAF(XMIN, 'SCALE' , XMAX, YMIN, 'SCALE' , YMAX) 
CALL CURVE(X,Y,200,0) 

CALL ENDPL(0) 

CALL DONEPL 


RETURN 
END 


SUBROUTINE RANGE(Y, YMIN, YMAX,N) ; 


Q aag 


DIMENSION Y(N) 
YMIN=1, E20 
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YMAX=- 1, E20 
DO 10 I=1,N 
IF(Y(1I).GT. YMAX) YMAX= Y(T) 
IF(Y(1). LT. YMIN) YMIN= Y(I) 
10 CONTINUE 


RETURN 
END 


SUBROUTINE RANGE1(Y,YMIN, YMAX,N) 


DIMENSION Y(N) 
YMIN=1. E20 
YMAX=-1. E20 
DO 20 I=1,N 
IF(Y(1).GT. YMAX) YMAX= Y(I) 
IF(Y(1). LT. YMIN) YMIN= Y¥(I) 
20 CONTINUE 


a aan 


RETURN 
END 
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APPENDIX B. B MAX - RISE TIME PROGRAM 
PROGRAM CERENKV 


Petetede dete te tete ele he Tee beTeHe Teter PAIR RRR RI RRA AAA AIA ARIANA AIAN IRER 


* 
* 
* 
* 
* 


~ 


ie 


+ 


THIS PROGRAM IS LOCATED ON THE NPS MAIN FRAME COMPUTER AND IS 
WRITTEN USING WATFOR 77. AFTER ENTERING CERTAIN PARAMETERS, THE 
SHAPE OF THE CERENKOV PULSE IS GENERATED. THE MAXIMUM VALUE OF 
THE PULSE MAGNITUDE IS DETERMINED AND STORED WITH ITS CORRES- 
PONDIND RISE TIME VALUE. THE RISE TIME IS INCREMENTED AND THE 
ENTIRE PROCESS IS REPEATED. THESE VALUES ARE THEN READ TO AN 
OUTPUT FILE FOR USE IN A PLOTTING ROUTINE. 


te tet + 


Letetetedeleketetete dete dede FededeLedekedzedne PALA IAI IRMA AARAALAS AAAS SAAD AI ADIDAS ASIA 


REAL N, U1, U2, BETA, CO, ROE, Al, A2, BPRME, YMAX 

REAL RISETM, D, DD, El, E2, Z1I, Z1F, Z21, Z2F 

REAL Bi, B2, W1, W2, YY, XX, T1, T2, WPI, WPF 

REAL ZPMIN, SSQ, BPSQ, TINC, IAMP, ENER, ENERJ 

REAL S, Z, C, MO, V, NEG, TPRME, B, DRISE, BMAX 

INTEGER I, J 

DIMENSION TPRME(200), B(200), DRISE(300), BMAX(300) 

DATA TPRME/200*0.0/, B/200*0.0/, DRISE/300*0.0/, BMAX/300*0. 0/ 


OPEN (UNIT=7 ,FILE='PULPLT DATA A') 
OPEN (UNIT=4,FILE='TERMINAL' ) 
OPEN (UNIT=7 ,FILE='CPLOT1 DATA A‘) 


PERRIER RTT NTT RAIL NNR IRIAN IRI IK IN IR IR IDR INI II RII III NII IKI IIE. 


X 
cy 


* 
* 


N IS THE INDEX OF REFRACTION. IAMP IS THE ACCELERATOR CURRENT * 
IN AMPS. ENER IS THE ELECTRON ENERGY IN MEV. RISETM IS THE * 
PULSE RISETIME IN SECONDS. S IS THE OBSERVER RADIAL DISTANCE IN * 
METERS. Z IS THE OBSERVER AXIAL DISTANCE IN METERS. CO IS THE * 
VACUUM SPEED OF LIGHT IN METERS PER SECONDS. MO IS THE ELECTRON * 
MASS IN KILOGRAMS. ROE IS THE LINEAR CHARGE DENSITY IN COULOMBS * 
PER METER. * 


Fetvrerederetereaededededetededetended de tetedevedetedeKe tN FRM RINT NTR LTTE 


N = 1.000293 
TAMP = 510.0 
ENER = 30.0 
RISETM = 9. 9E-12 
S = 10.0 

Z= 0.0 


CO = 2.99792458E08 
MO = 9.109534E-31 

C = CO/N 

ENERJ = ENER*1. 6021892E~-13 

V = CO*SQRT(1. O-( (MO%*2. 0) *( COMA, 9) / CENERI**2. 0))) 
IF (V.LE.C) GO TO 80 

BETA = V/CO 

ROE = IAMP/V 
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BPRME = N*BETA 


HARRI ATA TIT ITN NNT RTT RII RTI RIN IKARIA IIR III RR 


* SETTING Ul ARBITRARILY TO ZERO FIXES U2. THIS IS TO SIMPLIFY * 
* CALCULATIONS. * 
eA IIIA HITTITE ITH ITNT RI RNIN TH TR TTI IN TIT TT T TTH ET TR EH TET TE TIE 


U1 = 0.0 


Hemekedesedere dete Week Te te HK KARAT IMIR TT TRITON OO TCT TE TET TE TE TE TET ETE 


* WITH 200 POINTS, TINC WAS SELECTED SO AS FOR THE PULSE MAXIMUM * 
* 


* TO OCCUR AT APPROXIMATELY 20% OF THE GRAPH LENGHT. 
HAKATA AAR TRAIN TRT RIT TTIIR TNT RTT RITTER TT TER ET TTI TT TN THE 


SSQ = S**2,0 

BPSQ = BPRME**2. 0 
D = Z*BPSQ 

DD = 1 = BPSQ 


NEG = SSQ*DD 


Seteveretedededevedededevedcdelesededetedese dese sere vedere vee TTT TET TTR NCTE ETE TE ETE TET IE, 


* ZPMIN IS THE Z VALUE OF THE MINIMUM OF THE U CURVE. * 
aviedevereveveseded ede easededevededesedeasesevedeaevededenetevedbirdedk kevesk rede devesesedevevedesede ved ae lesese dete severe we 


ZPMIN = Z - (S/(SQRT(BPSQ - 1.0))) 


deskeredeakededekdeinbaiededededdedeekiedeiicleibicileciciekiedcdedeiendededbdededeebievedcaededeveteacdete sete neve 


* T1 AND T2 ARE THE TIMES THAT THE U CURVE MINIMUM IS TANGENT TO * 
* THE Ul AND U2 LEVELS RESPECTIVELY. SINCE T2 AND U2 ARE FUNCT- * 
* 


* IONS OF RISETM, THEY ARE INSIDE THE DO LOOP. 
devededsdedevededededesededededededevedede rede dete re dete reser TENNENT THEN TE TT TNT TE eT ITE Te ese TENCE 


T1 = (ZPMIN ~ U1)/V + (SQRT(SSQ + (Z - ZPMIN)**2.0))/C 


DO 60 J = 1,300 
RISETM = RISETM + 1, 0E-13 
TINC = RISETM/40. 0 
DRISE(J) = RISETM 
U2 = -V*RISETM 
T2 = Tl + RISETM 


DO 40 I = 1,200 
IF (I.EQ.1) THEN 
TPRMEC(I) = Ti 


B(I) = 0.0 
ELSE 

TPRME(I) = T1 + (REAL(1I)*TINC) 
END IF 


FetededetedetedetedeledeledeVeleledlese Fee L LAN TILA SSA SS SANE LSS SAAT A ASN ALN LARA AA RIT 


* THIS IF STATEMENT HANDLES THE CHANGEOVER IN LIMITS OF INTEGRAT- * 
*% ION FROM ONE INTERVAL TO TWO INTERVALS. % 


Tacdedevededetedededevedededede sede de Tete Sede dene reK TL RTI TT TET TE TE TENE LED INE TEL TLD SEDAN INDIA IAN 


IF (TPRME(I).GT.T2) GO TO 10 
Al = Ul + V*TPRMEC1) 
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El = (Z = Al)**2.0 


Feredededetede Heke Never Ne HITTIN ETI RTI TIRE RII NITION TOE 


* THIS IF STATEMENT HANDLES THE CASE OF HAVING A NEGATIVE SQUARE * 


* ROOT IN REAL SPACE WHEN DETERMINING THE LIMITS OF INTEGRATION. * 
HAART NT IRI NIT RTT IAN IRI TTI ITT TT RINT I TTT TEN TEN TET ITER TET TET CT HIE 


10 


IF (E1. LE. ABS(NEG)) THEN 
B(I) = 0.0 
GO TO 


(A1-D) + SQRT(BPSQ*(E1+NEG) ))/DD 
(A1-D) - SQRT(BPSQ*(E1+NEG) ))/DD 


WPI = Z - ZiT 
= Z- Z21F 
WPI/S 
WPF/S 
ATAN(W1) 
ATAN(W2) 
BCI) = -ROE*N*( BETA**2. 0)*(YY-XX)*1. 0E04/U2 
GO TO 20 


< 
< 
iu i 


CONTINUE 
Al = Ul + V*TPRME(I) 


A2 = U2 + V*TPRME(T) 
El = (@ - Al)**2,0 
E2 = (@ ~ A2)**2,0 


IF (E1. LE. ABS(NEG)) THEN 
B(I) = 0.0 
GO TO 30 


1-D) + SQRT(BPSQ*(E1+NEG) ))/DD 
1-D) - SQRT(BPSQ*(E1+NEG)))/DD 


aN 
tant 
> > 


B(I) = 0.0 


= ((A2-D) + SQRT(BPSQ*(E2+NEG) ))/DD 
= ((A2-D) - SQRTCBPSQ*(E2+NEG) )) /DD 


Z- Zul 

Z = Z21 

WPI/S 

WPF/S 

ATAN(W1) 

ATAN(W2) 

~ROE*N*(BETA**2, O)*(YY - XX)%1, OB04/U2 


< 
— 
nuenunu 


WPI = Z - 22F 
WPF = Z ~ ZiF 
Wl = WPI/S 
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W2 = WPF/S 

YY = ATAN(W1) 

XX = ATAN(W2) 

B2 = -ROE*N(BETA*#2. 0)*(YY - XX)*1. 0£04/U2 
B(I) = Bl + B2 


C 
20 CONTINUE 
C 
30 CONTINUE 
C 
40 CONTINUE 
C 
YMAX = +1. 0E20 
C 
Co Setibbeeddninhindieiitciset titi tint ePIC I TIT TI TH TI RTI IIIT IIIT 
CG %* THIS LOOP DETERMINES THE MAXIMUM VALUE OF THE CERENKOV PULSE. * 
Co Seitseinvineskininiebnibiiccicidcaiacindicdt aire ase r eee IIT TT RT ITI NII 
C 
DO 50 I = 1,200 
IF (B(1I).GT.YMAX) YMAX = B(TI) 
50 CONTINUE 
C 


BMAX(J) = YMAX 
60 CONTINUE 


DO 70 I = 1,300 
WRITE(7,1000) DRISE( 1) ,BMAX(I) 
1000  FORNAT(G15. 7,G15. 7) 
70 CONTINUE 
C 
GO TO 90 
80 WRITE(4,1010) 
1010 FORMAT(/,'V IS LESS THAN C, THEREFORE THERE IS NO CERENKOV RADIATI 
$ON FOR THIS ELECTRON’ ,4X,' ENERGY’) 
90 CONTINUE 


STOP 
END 


PROGRAM CPLOT1 
PROGRAM LOGPLT 


a aga aang 


REAL X(300),¥(300) ,XMIN, XMAX, YMIN, YMAX 

DOUBLE PRECISION X(300) ,¥(300) ,XMIN,XMAX, YMIN, YMAX 
INTEGER N 

CHARACTER*20 TITLES 


TITLES = 'B MAX VS RISETIMES' 
N = 300 


READ(7,8000) ((X(I),Y(I)}),I=1,N) 
8000 FORMAT( 2615. 7) 


C CALL TEK618 
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10 


20 


CALL SHERPA('P7 he ce 

CALL SHERPA('P8 , vip? s 

CALL PAGE(6,7.5) 

CALL AREAZD(4. 5,5.5) 

CALL XNAME(' RISE TIME (SECONDS) $' »100) 
CALL YNAME('CERENKOV PULSE (GAUSS)5' ,100) 
CALL HEADIN(TITLES,100,1.5,1) 

CALL CROSS 

CALL RANGE(X,XMIN, XMAX,N) 

CALL RANGE1(Y, YMIN, YMAX , N) 

CALL GRAF(XMIN, ‘SCALE! SXMAX sYMIN, 'SCALE' , YMAX) 
CALL LOGLOG(XMIN,4. 5 SYMIN, 5.5) 

CALL CURVE(X,Y,300,0) 

CALL ENDPL(0) 

CALL DONEPL 


RETURN 
END , 


SUBROUTINE RANGECY, YMIN, YMAX,N) 


DIMENSION Y(N) 
YMIN=1. E20 
YMAX=~1. £20 
DO 10 I=1,N 
IF(Y(1).GT. YMAX) YMAX= Y(1) 
IF(Y(I). LT. YMIN) YMIN= Y(1) 
CONTINUE 


RETURN 
END 


SUBROUTINE RANGE1(Y, YMIN, YMAX,N) 


DIMENSION Y(N) 
YMIN=1. E20 
YMAX=-1. £20 
DO 20 I=1,N 
IF(Y(1).GT. YMAX) YMAX= Y(T) 
IF(Y(1). LT. YMIN) YMIN= Y(I) 
CONTINUE 


RETURN 
END 
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APPENDIX C. BMAX -S PROGRAM 
PROGRAM BVS 


PRAIA ARAN A ATARI ARITA IAAI TAA TIAA NIA DIARIO IISA IAA AIA AI IR 
* THIS PROGRAM IS LOCATED ON THE NPS MAIN FRAME COMPUTER AND IS * 
WRITTEN USING WATFOR 77. AFTER ENTERING CERTAIN PARAMETERS, THE * 
SHAPE OF THE CERENKOV PULSE IS GENERATED. THE MAXIMUM OF THE * 
PULSE MAGNITUDE IS DETERMINED AND STORED WITH ITS CORRESPONDING * 
S VALUE. THE S VALUE IS THEN INCREMENTED AND THE ENTIRE PROCESS * 
IS REPEATED. THESE VALUES ARE READ TO AN OUTPUT FILE FOR LATER i 


USE IN A PLOTTING ROUTINE. 
HIRATA RAAT ATTA IAAI AI IARI TAIRA AAS IIIA RL LIA IIIA IIIT III IIH 


+t 4 t + 


REAL N, U1, U2, BETA, CO, ROE, Al, A2, BPRME, YMAX 

REAL RISETM, D, DD, E1, E2, Z1I, Z1F, 221, Z2F 

REAL Bl, B2, W1, W2, YY, XX, T1, T2, WPI, WPF 

REAL ZPMIN, SSQ, BPSQ, TINC, IAMP, ENER, ENERJ 

REAL S, Z, C, MO, V, NEG, TPRME, B, DS, BMAX 

INTEGER I, J 

DIMENSION TPRME( 200), B(200), DS(300), BMAX(300) 

DATA TPRME/200*0.0/, B/200%*0.0/, DS/300*0.0/, BMAX/300*0, 0/ 


OPEN (UNIT=4,FILE='TERMINAL' ) 
OPEN (UNIT=7,FILE='BVSPLT DATA A‘) 


devevedededededevedededededetedetedetete terete selene Ne Teh Tete HLT L LTE RVR TRUDE RSD SELLS LITLE L DSIRE 
*%* N IS THE INDEX OF REFRACTION. IAMP IS THE ACCELERATOR CURRENT * 
%* IN AMPS. ENER IS THE ELECTRON ENERGY IN MEV. RISETM IS THE ve 
%* PULSE RISETIME IN SECONDS. S IS THE OBSERVER RADIAL DISTANCE IN * 
* METERS. Z IS THE OBSERVER AXIAL DISTANCE IN METERS. CO IS THE * 
%* VACUUM SPEED OF LIGHT IN METERS PER SECONDS. MO IS THE ELECTRON * 
%* MASS IN KILOGRAMS. ROE IS THE LINEAR CHARGE DENSITY IN COULOMBS * 
* PER METER. ve 
Sevedevevedededededetedede Teste Pere Te eT TRIN NTT NTE ELIE ST RAEL ER IREN II NSIS IN INVI NAG 


N = 1, 000293 
TAMP = 510.0 
ENER = 30.0 
RISETM = 10. 0E-12 
5 = 0.0 
Z = 0.0 


CO = 2.99792458E08 
MO = 9. 109534E-31 

C = CO/N 

ENERJ = ENER*1, 6021692E-13 

V = CO*SQRT( 1. 0-( (MO**2. 0)*(CO**4, 0) /(ENERI**2. 0))) 
IF (V.LE.C) GO TO 80 

BETA = V/CO 

ROE = IAMP/V 

BPRME = N*BETA 
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aagaaaga aangnana aagaaAaaAan 


aQaaana 


ana 


MAHA IRN RTT RIT TTI ITT TTR ITT RIC AR ITI EON TOE 


* SETTING Ul ARBITRARILY TO ZERO FIXES U2. THIS IS TO SIMPLIFY * 
* CALCULATIONS. = 
frininiiciininciinictinindictet tear ree eae AIT ITT TTT IIT TTT I TR TTT TE 


U1 
U2 


0.0 
-V*RISETM 


PRR IRRITATED REPRE RRR I ARI IIINT A N IRI ARIRAATOIAIR 


* WITH 200 POINTS, TINC WAS SELECTED SO AS FOR THE PULSE MAXIMUM 


%* TO OCCUR AT APPROXIMATELY 20% OF THE GRAPH LENGHT. 
Pa Loire arg eases ach aa noe ein ree Re ETI 


TINC = RISETM/40. 0 
BPSQ = BPRME**2, 0 
D = Z*BPSQ 

DD = 1 - BPSQ 


DO 60 J = 1,300 
$=8§ + 0,01 
DS(J) = 

SSQ = S**2.0 

NEG = SSQ*DD 


aedetedetedederedesetedede Reese vesete eee RATNER TT TTT TNR TT TONITE NIT NC TN TENT TOT TENET TET ETON TOTO 


* ZPMIN IS THE Z VALUE OF THE MINIMUM OF THE U CURVE. * 
HAART NTR TEIN TEN TEFEN TET TEIN TET NTN TITRE INI TTR TERETE TEE TE FEN TEE TE ETE 


ZPMIN = 2 - (S/(SQRT(BPSQ - 1.0))) 


Hietevesedededesetededededek deter sete Ne sewe Fe TeTeR eM TeV Te MN INTC TTR INTENT TTI TAT TCTE NT TN IRI, 


* T1 AND T2 ARE THE TIMES THAT THE U CURVE MINIMUM IS TANGENT TO * 
* THE Ul AND U2 LEVELS RESPECTIVELY. * 
Sedededededetededededetedetetele ede leseleledele Lele TR NTL TALL NNSA INNA ARIEL IIASA INNIS IIIS 


Tl 
T2 


(ZPMIN - U1)/V + (SQRT(SSQ + (Z - ZPMIN)**2.0))/C 
Tl + RISETM 


DO 40 I = 1,200 
IF (I.EQ.1) THEN 
TPRME(I) = Tl 


B(I) = 0.0 
ELSE 

TPRME(I) = T1 + (REAL(1I)*TINC) 
END IF 


RATA TINIAN NT TRI TRELI TTL TRTRTRNCRIR TR IT TR INITTRTER ERE RET 


% THIS IF STATEMENT HANDLES THE CHANGEOVER IN LIMITS OF INTEGRAT- * 
* 


* ION FROM ONE INTERVAL TO TWO INTERVALS. 
Peledededede dele deLededePeleheFeLedede lV LVL TRL ALAA LTA IIS SES SASL NAA AEA AILS SISA AAAI 


IF (TPRME(1I).GT.T2) GO TO 10 
Al = Ul + V*TPRME(I) 
= (Z + Al)**2.0 


devedededededede dete de Lelede le ede Le BeLede TeV ERL ERENT 
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* THIS IF STATEMENT HANDLES THE CASE OF HAVING A NEGATIVE SQUARE * 


C 

C * ROOT IN REAL SPACE WHEN DETERMINING THE LIMITS OF INTEGRATION. * 
: Co ininintveideisdirineiieinedsbicaeiesteacaeaeser i neven sen rene rere ae Wie TNT THT HT TET TIT TIE 

C 


IF (E1. LE. ABS(NEG)) THEN 
B(I) = 0.0 
GO TO 30 


(Al-D) + SQRT(BPSQ*(E1+NEG)))/DD 
(A1-D) = SQRT(BPSQ*(E1+NEG) ))/DD 


WPI = Z - Z1I 

WPF = 

W1 = WPI/S 

W2 = WPF/S 

YY = ATAN(W1) 

XX = ATAN(W2) 

BCI) = -ROE*N*( BETA**2, 0)*(YY-XX)*1. 0E04/U2 
GO TO 20 





10 CONTINUE 
Al = Ul + V*TPRME(T) 
= U2 + V*TPRME(I) 
El = (Z - Al)**2.0 
E2 = (Z - A2)**2.0 
IF (E1. LE. ABS(NEG)) THEN 
B(I) = 0.0 


Z1I = ((A1-D) + SQRT(BPSQ*(E1+NEG) ))/DD 
= ((AL-D) + SQRT(BPSQ*(E1+NEG) ))/DD 


IF (E2. LE. ABS(NEG)) THEN 
B(I) = 0,0 
GO TO 30 


(CA2-D) + SQRT(BPSQ*(E2+NEG) )) /DD 
(CA2-D) - SQRT(BPSQ*(E2+NEG) ))/DD 


Z- 211 

B= °Z21 

WPI/S 

WPF/S 

ATAN(W1) 

ATAN(W2) 

-ROE*N*( BETA**2. 0)*(YY - XX)"1. 0£04/U2 


erg 
wuuwu nh 


WPI = Z - 22F 
=2- Z1F 
WPI/S 
WPE/S 
ATAN(W1) 
ATAN(W2) 


? 

re 

< 
Huu lt 














B2 = -ROE*N*(BETA**2,0)*(YY = XX)*1, OE04/U2 
B(I) = Bl + B2 


C 
20 CONTINUE 
C 
30 CONTINUE 
C 
40 CONTINUE 
C 
YMAX = +1. 0E20 
DO 50 I = 1,200 
IF (B(I).GT. YMAX) YMAX = B(I) 
50 CONTINUE 
BMAX(J) = YMAX 
60 CONTINUE 
C 


DO 70 I = 1,300 
WRITE(7,1000) DS(I) ,BMAX(I) 
1000 FORMAT(G15. 7,G15. 7) 
70 CONTINUE 


GO TO 90 

80 WRITE(4,1010) 

1010 FORMAT(/,'V IS LESS THAN C, THEREFORE THERE IS NO CERENKOV RADIATI 

SON FOR THIS ELECTRON' ,4X,' ENERGY’) 

90 CONTINUE 

Cc 

STOP 
END 


PROGRAM BVSPLT 


REAL X(300) ,Y(300) ,XMIN,XMAX, YMIN, YMAX 

DOUBLE PRECISION X(300),¥(300) ,XMIN,XMAX, YMIN , YMAX 
INTEGER N 

CHARACTER*25 TITLES 


a Q aaa 


TITLES = 'B MAX VS RADIAL DISTANCES' 
N = 300 


READ(7,8000) ((X(1),Y(1I)),I=1,N) 
8000 FORMAT(2G15. 7) 


c CALL TEK618 
CALL SHERPA('P9 ",'B',3) 
CALL PAGE(6,7.5) 
CALL AREA2D(4.5,5.5) 
CALL XNAME('RADIAL DISTANCE (METERS)$' ,100) 
CALL YNAME('CERENKOV PULSE (GAUSS)$' ,100) 
CALL HEADIN(TITLES,100,1.5,1) 
CALL CROSS 
CALL RANGE(X,XMIN, XMAX,N) 
CALL RANGE1(Y,YMIN, YMAX,N) 
CALL GRAF(XMIN, ‘SCALE’ ,XMAX,YMIN, ‘SCALE’ , YMAX) 
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CALL CURVE(X,Y,300,0) 
CALL ENDPL(0) 
CALL DONEPL 


RETURN 
END 


SUBROUTINE RANGE(Y,YMIN, YMAX,N) 


DIMENSION Y(N) 
YMIN=1, E20 
YMAX=-1. E20 
DO 10 I=1,N 
IFCY(I). GT. YMAX) YMAX= Y(1I) 
IF(Y(1). LT. YMIN) YMIN= Y(I) 
CONTINUE 


RETURN 
END 


SUBROUTINE RANGE1(Y, YMIN, YMAX,N) 


DIMENSION Y(N) 
YMIN=1, E20 
YMAX=-1, E20 
DO 20 I=1,N 
IF(Y(1I). GT. YMAX) YMAX= Y(1I) 
IF(Y(I). LT. YMIN) YMIN= Y(1I) 
CONTINUE 


RETURN 
END 
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APPENDIX D. S-T-B PROGRAM 
PROGRAM STB 


WRN AAR AANA AHT ARI ARIA AH AIHA IARI IRA IIA IAAI TAIN IIR 
* THIS PROGRAM IS LOCATED ON THE NPS MAIN FRAME COMPUTER AND IS * 
* WRITTEN USING WATFOR 77. AFTER ENTERING CERTAIN PARAMETERS, THE * 
* SHAPE OF THE CERENKOV PULSE IS GENERATED. THE S VALUE IS INCRE- * 

¥% MENTED AND THE PROCESS IS REPEATED. THESE VALUES ARE STORED FOR 7 


* LATER USE IN A 3-D PLOTTING ROUTINE. 
Se EE eee EL 


REAL N, U1, U2, BETA, CO, ROE, Al, A2, BPRME 
REAL RISETM, D, DD, E1, E2, Z1I, ZF, 221, 22F 
REAL Bl, B2, Wl, W2, YY, XX, Tl, T2, WPI, WPF 
REAL ZPMIN, SSQ, BPSQ, TINC, IAMP, ENER, ENERJ 
REAL S, Z, C, MO, V, NEG, TPRME, B 

REAL BF, SRD, TIME 

INTEGER I, J, K, L, M 

DIMENSION TPRME(200), B(200) 

DIMENSION BF(2500), SRD(2500), TIME(2500) 

DATA TPRME/200%0,0/, B/200*0.0/, M/2500/ 

DATA BF/2500%*0.0/, SRD/2500%0.0/, TIME/2500*0. 0/ 


OPEN (UNIT=4,FILE='TERMINAL' ) 
OPEN (UNIT=7,FILE='S DATA B') 
OPEN (UNIT=8 ,FILE='T DATA B') 
OPEN (UNIT=9,FILE='B DATA B') 


FeWeHe Hee WeNeWeweKe Herne Wwe dew ye ne ese Fe He ve dese de KVM MIN LAIR IN, BINT ITE TENE TEER ATEN IEE NITE IE 
* N IS THE INDEX OF REFRACTION. IAMP IS THE ACCELERATOR CURRENT * 
%* IN AMPS. ENER IS THE &£LECTRON ENERGY IN MEV. RISETM IS THE * 
* PULSE RISETIME IN SECONDS. S IS THE OBSERVER RADIAL DISTANCE IN * 
* METERS. 2 IS THE OBSERVER AXIAL DISTANCE IN METERS. CO IS THE * 
% VACUUM SPEED OF LIGHT IN METERS PER SECONDS. MO IS THE ELECTRON * 
MASS IN KILOGRAMS. ROE IS THE LINEAR CHARGE DENSITY IN COULOMBS ‘i 


PER METER, 
cirri DETTE Te TT eI ie eee RET 


+ % 


N = 1.000293 

= 510.0 
ENER = 30.0 
RISETM = 10. 0OE-12 

0. 

0 


= 2.99792458E08 

= 9,109534E-31 

CO/N 

= ENER*1. 6021892E-13 

CO*SQRT( 1. O~( (MO**2, 0) %*( CO**4, 0) /(ENERJ**2.0))) 
(V.LE.C) GO TO 90 

ETA = V/CO 


WHF aQ NM 
A ao 
1] Ss a || 
7 
aq 


ay 
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aqgaangagIga aQaaaqaan 


ROE = JAMP/V 
BPRME = N*BETA 


HANH HITHER ANAT NTI IRI NTT TI THT IR ITN TIT IT TER TTT IIE 
* SETTING U1 ARBITRARILY TO ZERO FIXES U2. THIS IS TO SIMPLIFY * 
* 


* CALCULATIONS. 
Aainrinie inane teiten venir ire eren reer IER TTI RI TOT ET ER TTT A 


Ul 


0.0 
U2 = -V*RISETM 


eR AANA ANIA TAT KART RIT RR RTI IRIN LIRIK CRITI RRR IIR. 


%* WITH 200 POINTS, TINC WAS SELECTED SO AS FOR THE PULSE MAXIMUM * 
%* 


* TO OCCUR AT APPROXIMATELY 20% OF THE GRAPH LENGHT. 
Petederetesedetedede dele vested Tee te tere sede Bese KVM TTT TTT ICR TOT TE TE TENE TENE TE TE TORE TOE TE NTE TENE TE TENE TENE 


TING = xXISETM/10. 0 


BPSQ = BPRME**2. 0 
D = Z2*BPSQ 
DD = 1 - BPSQ 


DO 50 K = 1,50 
IF (K.EQ.1) THEN 
J=0 


ELSE 

J = 50*(K-1) 
END IF 
§=8§ +0,.5 
SSQ = S20 
NEG = SSQ*DD 


Feverlededederevedevedededededetededte dese deve Pevere dese eek TEL LTT TELE TEP NLL LTD NLS LER IAIN LER VIERA 
%* ZPMIN IS THE Z VALUE OF THE MINIMUM OF THE U CURVE. * 


Pebeledededededese deste de ele Leelee LT LTTE LLANES AISLE INAS AAAI NAINA ADRE 


ZPMIN = Z@ - (S/(SQRT(BPSQ - 1.0))) 


etededededevedesedededededetede tele deeded cevele EINSTEIN NSD NENA SD IID SESS BINS SSA SNAILS 
%* Ti AND T2 ARE THE TIMES THAT THE U CURVE MINIMUM IS TANGENT TO * 
% 


* THE U1 AND U2 LEVELS RESPECTIVELY. 
Terederertodededevede dete tedesetese te dee vee We TeV RITE TNH TENTION TEE TENE TELE TE TE TENET NE TENE TENE TE TENE TEN TENE NEE 


ZPMIN = U1)/V + (SQRT(SSQ + (Z - ZPMIN)**2.0))/C 
1 + RISETM 


i an 


DO 40 I = 1,50 
IF (I.EQ.1) THEN 
TPRME(I) = Tl 


BCI) = 0. 


y 
3 
is? | 
oN 
= 
ww 
tt © il 


T1 + (REAL(1)*TINC) 


SRD(L) = § 
TIME(L) = TPRME(1) 














qQaaac 


AaAaINIIaAo 


HeIeiekelek tei rH sew dw—enwwen wen Wsrez— ne dwegzcw7ez KaARIA IAAT IRI IRR TAI RIKI TTT ITE TET 


%* THIS IF STATEMENT HANDLES THE CHANGEOVER IN I:IMITS OF INTEGRAT- : 


* TON FROM ONE INTERVAL TO TWO INTERVALS. 
HIT MINIR ANAT HRN AT HIRT IHN I HTT HTT IN IITA III IIR RITTER ET 


IF (TPRME(1I).GT.T2) GO TO 10 
Al = U1 + V*TPRME(T) 
El = (Z = A1)**2.0 


HAAR RAR ITN IR IARI RNAI ARIANA ARIAS EPRI APPR RIT TINTTTETE 


* THIS IF STATEMENT HANDLES THE CASE OF HAVING A NEGATIVE SQUARE * 


* ROOT IN REAL SPACE WHEN DETERMINING THE LIMITS OF INTEGRATION. * 
dedesereredeteedtee dese AIR HII RT RT IR TNR TR TTR IR IRIN TORT TET TET TORT TCT TI 


IF (E1. LE. ABS(NEG)) THEN 


B(I) = 0.0 
GO TO 30 
RLSE 


Z1I = ((A1-D) + SQRT(BPSQ*(E1+NEG)))/DD 
Z21F = ((Al-D) «+ SQRT(BPSQ*(E1+NEG) ))/DD 
ENDIF 


WPI = 4 - ZilI 

WPF = Z ~ ZiF 

Wi = WPI/S 

W2 = WPF/S 

YY = ATAN(W1) 

XX = ATAN(WZ) 

BCI) = -ROE*N*(BETA**2. 0)**(YY-XX)*1. OE04/U2 
GO TO 20 


10 CONTINUE 

Al = U1 + V*TPRME(I) 

A2 = U2 + V*TPRME(1) 

El = (Z ~ Al)**2.0 

E2 = (Z = A2)**2.0 

IF (E1. LE. ABS(NEG)) THEN 
B(I) = 0.0 
GO TO 30 


1-D) + SQRT(BPSQ*(E1+NEG) ))/DD 
1-D) - SQRT(BPSQ*(E1+NEG) )) /DD 


N 
_ 
eRe 
oi 
NIN 
oNoN 
> > 


IF (£2. LE. ABS(NEG)) THEN 

BCI) = 0.0 

GO TO 30 
ELSE 

Z2I = ((A2-D) + SQRT( BPSQ*(E2+NEG) ) )/DD 
Z2F = ((A2-D) - SQRT(BPSQ*(E2+NEG) ))/DD 
ENDIF 


WPI = 2 - 211 ; 
WPF = Z = 221 

W1 = WPI/S 

W2 = WPF/S 
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20 
30 


40 
50 


80 
C 


90 
1010 


100 


qQ 


Qa aaa 


YY = ATAN(W1) 

XX = ATAN(W2) 

Bl = =ROE*N*( BETA**2, 0)*(YY - XX)*1. 0E04/U2 
WPI = Z + Z2F 

WPF = Z - Z1F 

W1 = WPI/S 

W2 = WPF/S 

YY = ATAN(W1) 

XX = ATAN(W2) 


B2 = -ROE*N*(BETA**2,0)*(YY - XX)*1. 0E04/U2 
B(I) = Bl + B2 


CONTINUE 
CONTINUE 


BF(L) = B(T) 
CONTINUE 


CONTINUE 


DO 60 I = 1,M 
WRITE(7,1000) SRD(I) 
FORMAT(G15. 7) 

CONTINUE 


DO 70 1 = 1,M $ 
WRITE(8,1000) TIME(T) 
CONTINUE 


DO 80 I = 1,M 
WRITE(9,1000) BF(1) 
CONTINUE 


GO TO 100 
WRITE(4, 1010) 
FORMAT(/,'V IS LESS THAN C, THEREFORE THERE IS NO CERENKOV RADIATI 


SON FOR THIS ELECTRON' ,4X,' ENERGY') 


CONTINUE 


STOP 
END 


PROGRAM THREED 

REAL 8(2500) ,T( 2500) ,B( 2500) ,X3MAX, Y3MAX , Z3MAX , ZMAT( 2500) 
INTEGER I,N 

CHARACTER*35 TITLE$ 


TITLES = 'B-FIELD VS TIME VS RADIAL DISTANCES' 
N = 2500 


READ(7,8000) (SCI), I=1,N) 
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Q aaa 


8000 FORMAT(G15. 7) 
READ(8,8000) (T(I),I=1,N) 
READ(9,8000) (BCI) ,I=1,N) 


CALL RANGE(S ,X3MAX,N) 
CALL RANGE(T, Y3MAX,N) 
CALL RANGE(B ,Z3MAX,N) 


CALL TEK618 

CALL SHERPA('P10 "5 'B';3) 

CALL PAGE(7.5,9. 0) 

CALL AREA2D(5.5,7.0) 

CALL VOLM3D(1,1,1) 

CALL X3NAME('DISTANCE (METERS)$"‘, 100) 
CALL Y3NAME('TIME (SECONDS)$' , 100) 
CALL Z3NAME('MAGNITUDE (GAUSS)$',100) 
CALL HEADIN(TITLE$,100,1.0,1) 

CALL VUABS(1.0,-0.5,0.5) 

CALL GRAF3D(0. 0, 'SCALE' ,X3MAX,0.0,'SCALE',Y3MAX,0. 0, 'SCALE' , 23MAX) 
CALL BGNMAT(50,50) 

CALL GETMAT(S,T,B,2500,0) 

CALL ENDMAT( ZMAT, 0) 

CALL SURMAT( ZMAT,1,50,5,50,0) 

CALL ENDPL(0) 

“ALL DONEPL 


STOP 
END 


SUBROUTINE RANGE(Y, YMAX,N) 


DIMENSION Y(N) 
YMAX = -1.0E20 
DO 20 I = 1,N 
IF (Y(1).GT. YMAX) YMAX = Y(T) 
20 CONTINUE 


RETURN 
END 
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